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4. 巴 基 


f X: 当前 ,基于 正 积温 的 度 日 模型 广泛 


因子 ,并 分 析 了 度 日 因子 的 时 空 变化 特 4 
别 为 2.83 mm + d^! - 4C 3.74 mm: d! - Co 
海拔 递增 率 分 别 为 0.003 7 mm + d'' .CC - 


斯 坦 气 象 局 ,伊斯兰 


应 用 于 冰川 消融 研究 
游 Sachen , Gharko , Barpu 冰川 2014—2016 年 的 物质 平衡 和 气温 实 
征 及 影响 因素 。 研 究 结 果 显 

73.91 mn: d^! 
m^,0.007 4 mm + d^! 
m ;对 于 同一 观测 点 而 言 , 度 日 因子 不 是 一 个 常数 ,会 随 着 时 间 的 变化 而 改变 


f& 44000) 


中 ,该 模型 的 核心 参数 是 度 日 因子 。 根 据 印 度 河 上 
测 资料 ,计算 得 到 消融 期 内 各 冰川 研究 区 的 度 日 
显示 :Sachen、Gharko、Barpu 冰川 度 日 因子 均值 分 
0C ;各 冰川 度 日 因子 皆 随 着 海拔 升 高 而 递增 ， 

Ci.m!0.0041mm-.d C7. 

,冰川 度 日 因子 随 着 年 际 变化 呈 增 加 


的 趋势 ; 度 日 因子 受 
LEA si E 
一 定 的 影响 ,面向 阳 坡 的 冰川 度 日 


BH 
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冰川 是 气候 变化 的 指示 剂 ”” ,冰川 变化 能 清 
晰 地 反映 区 域 乃 至 全 球 气候 变化 趋势 " 。 物 质 平 
衡 是 反应 冰川 响应 气候 变化 的 重要 指标 ,因而 估 
算 和 重建 冰川 物质 平衡 具有 重要 的 科学 意义 ,一 直 
受到 国内 外 众多 研究 者 的 关注 “””。 冰 川 物质 平衡 
是 冰川 积累 和 消融 的 代数 和 ,其 中 冰川 的 消融 变化 


磺 影 响 显著 , 度 日 因子 整体 上 随 着 表 磷 厚度 的 增加 而 递减 。 
促进 了 冰川 的 消融 , 表 磺 覆盖 区 冰川 度 日 因 
因子 随 海拔 递增 率 大 了 


A 


然而 表 研 厚度 低 于 2 cm 时 , 表 


子 大 于 裸露 区 冰川 ;冰川 朝 向 的 变化 对 度 日 因子 产生 了 
明 坡 。 


基于 气温 与 冰川 消融 的 良好 统计 关系 建立 的 经 验 统 
计 模 型 … 。 度 日 模型 的 输入 数据 是 气温 ,数据 获 
取 容 易 ,计算 简单 ,虽然 不 能 揭示 冰川 物质 平衡 变化 
的 物理 机 制 ,但 在 区 域 或 流域 尺度 上 ,通常 能 获得 与 
能 量 平衡 模型 类 似 的 模拟 效果 。 因 此 在 冰川 物质 平 
衡 变 化 ,冰川 融 水 估算 等 研究 中 , 度 日 模型 相对 于 能 


特征 更 能 反映 冰川 对 气候 变化 的 响应 程度 。 因 
此 ,全 面 开 展 冰川 物质 平衡 变化 尤其 是 ok ITNE 


化 研究 ET BA A Ok) Fak VER HS OR LT 
至 全 球 气候 系统 变化 特征 具有 重要 的 科学 价值 和 社 
会 经 济 意义 。 

国内 外 已 有 的 冰川 物质 平衡 研究 中 ,基于 能 量 
平衡 的 冰川 消融 模型 和 基于 气温 的 统计 模型 应 用 广 
泛 ""。 前 者 理论 基础 和 模型 结构 较 完 善 , 能 从 物理 
机 制 的 角度 有 效 地 揭示 冰川 消融 过 程 "”。 然 而 能 
量 平衡 模型 输入 参数 较 多 ,野外 观测 易 受 基础 条 件 


量 平衡 模型 应 用 更 为 广泛 。 
印度 河 发 源 于 青藏 高 原 西 部 西藏 阿里 地 区 的 冈 
底 斯 山 , 流 经 克什米尔 进入 巴基斯坦 ,于 巴基斯坦 卡 
拉 奇 人 阿拉 伯 海 。 其 上 游 流 域 地 跨 兴 都 库 什 一 喀 喇 
昆仑 一 喜马拉雅 (HKH) 三 大 著名 山脉 ,是 全 球 中 低 
纬度 冰川 发 育 ie ae ae 
ERA 11 413 条 冰川 ,覆盖 面积 达 15 061. 74 
p , BUR ENHISDEDITBUR REA T Ea 
ER Ee 
EREK, EA 20 世纪 中 期 以 来 , 兴 都 


限制 ,不 易 获取 足够 的 观测 数据 ,因此 能 量 平 衡 模型 
只 适用 于 观测 条 件 较 好 的 部 分 冰川 观测 区 '”) 。 后 
者 ,通常 被 称 为 “ 度 日 模型 "或 “ 度 日 因子 模型 ,是 
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或 微弱 前 进 , 称 之 为 “ 喀 喇 昆仑 异常 "2 。 这 说 明 
该 地 区 冰川 变化 异 于 全 球 其 他 地 区 ,学 术 界 对 于 
“ 喀 喇 昆仑 异常 "的 认识 存在 争议 ,有 待 未 来 更 多 深 
入 的 研究 。 因 此 ,选择 印度 河上 游 流域 典型 冰川 进 
行 水 文 气象 观测 ,开展 冰川 物质 平衡 变化 研究 对 于 
揭示 印度 河上 游 流域 冰川 变化 与 区 域 水 量 平衡 \ 流 
域 水 资源 分 配 乃至 区 域 经 济 发 展 具有 重要 的 意义 。 

本 文 所 研究 的 Sachen , Gharko , Barpu 冰川 都 积 
累 了 丰富 的 水 文 气象 观测 资料 。 本 文 从 各 典型 冰川 
的 观测 资料 人手, 分析 冰川 度 日 因子 的 空间 特征 及 
其 影响 因素 ,为 未 来 在 印度 河上 游 流 域 开展 冰川 融 
水 径流 人 研究 和 冰川 物质 平衡 估算 提供 度 日 因子 
参数 。 


1 研究 区 概况 与 数据 方法 


1.1 研究 区 概况 

本 研究 选取 印度 河上 游 流域 内 的 3 个 冰川 开展 
人 研究 工作 ,分 别 是 Barpu 冰川 .Gharko 冰川 和 Sachen 
冰川 (图 1)。 近 年 来 ,中 国 科学 院 青藏 高 原 研究 所 
(ITPCAS) 同 巴基斯坦 空间 与 大 气 层 大 气 研 究 委 员 
会 (Space and Upper Atmosphere Research Commis- 
sion) 巴基斯坦 气象 局 ( Pakistan Meteorological De- 
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partment ) 和 巴基斯坦 水 电 开 发 署 (Water and Power 
Development Authority ) 进行 了 紧密 的 合作 并 联合 开 
展 野 外 考察 ,对 印度 河上 游 流 域 选取 的 3 个 冰川 径 
流 和 物质 平衡 进行 了 观测 ,获取 了 各 监测 冰川 的 气 
象 水 文 数据 。 

Sachen 冰川 (35.329°N, 74. 752* E ) 位 于 西 喜 
马 拉 雅 山脉 北 坡 ,印度 河上 游 流 域 Astore 子 流域 境 
Wo 冰川 朝向 NE。 冰 川 面积 9. 89 km , v fifi E 
0.88 km ,全 长 8. 5 km, 平均 宽度 1.16 km ,海拔 
5 039 ~3 382 m, VK JI] ifti EE BEA T B) R 
盖 ,裸露 区 极 少 。 

Gharko 冰川 (36. 04?N, 74. 67* E ) 位 于 西 喀 喇 
昆仑 山脉 南 坡 ,印度 河上 游 流域 Gilgit 子 流域 境内 。 
冰川 朝向 SW。 冰 川 面积 30. 496 km?, 冰 储量 4. 03 
km? ,全 长 13.9 km ,海拔 6 906 ~2 953 m, 冰 川 自 海 
J£ 4 100 m 以 下 被 厚度 不 均 的 表 硬 覆盖 ,冰川 区 可 
分 为 高 海拔 裸露 冰川 区 和 表 厂 覆盖 冰川 区 。 

Barpu 冰川 (35.329°N, 74. 752?E ) 位 于 西 喀 喇 
昆仑 山脉 北 坡 , 印 度 河上 游 流域 Hunza 子 流 域 境 内 。 
该 冰川 是 由 Sumair 和 Miar 两 支 同 朝 北 向 的 分 支 冰 
川 在 海拔 3 300 m 附近 汇流 而 成 的 复式 山谷 
冰川 。 冰 川 朝向 NE 。 冰川 面积 88.248 km? , 冰 储 


36? 10' 0" N 


0 4500 9000 


18 000 m 
a | 

o 2014 年 消融 观测 点 
。 2015 年 消融 观测 点 
六” 自动 气象 站 

[—.] Barpu 冰 川 区 


74° 40' 0" E 


36*0' 0" N 


74° 50' 0" E "5*0 0" E 


0 1150 2300 4600m 


35? 20' 0" N 


2014 年 消融 观测 点 
2015 年 消融 观测 点 
2016 年 消融 观测 点 
自动 气象 站 
Sachenvk JI] [EX 


74° 44 0" E 74° 46’ 0" E 74° 48’ 0" E 


图 1 研究 区 位 置 及 花 杆 观测 点 分 布 示意 图 


Fig.1 Location of the study area and the distribution of observational sighting rods 


量 16.93 km ,全 长 24.5 km ,平均 宽度 3.6 km , 海 
1& 6 982 ~2 829 m, Sumair 冰川 以 及 冰川 汇流 区 上 
部 较 洁 净 , 表 矿 很 少 ,Miar 冰川 在 海拔 4 200 m 以 下 
至 冰川 汇流 点 区 域内 表 磺 较 多 ,4 200 m 以 上 部 分 
较 洁净 。 冰 川 汇流 区 在 海拔 3 100 m 以 下 直到 亲身 
末端 者 被 厚度 不 均 的 表 磺 覆盖 。 
1.2 数据 来 源 

3 个 冰川 的 消融 数据 来 自 2014 一 2016 年 的 夏 、 
秋季 野外 考察 结果 。 受 野外 条 件 的 限制 , 花 杆 主要 
布设 在 冰川 消融 区 ,冰川 消融 速率 较 快 , 强 年 多 数 花 
杆 完全 出 露 , 故 所 有 花 杆 的 布设 和 观测 均 在 观测 年 
内 完成 。 花 杆 观 测 内 容 包 括 各 布点 花 杆 处 的 海拔 、 
位 置 表 碘 厚 度 、 花 杆 出 圳 长度。 观测 期 内 ,所 有 花 
杆 无 积 雪 和 覆盖。 排除 观测 期 内 倒塌 或 折断 的 花 杆 ， 
2014—2016 年 观测 期 间 内 ,Sachen 冰川 累计 收获 了 
29 根 有 效 花 杆 , Gharko 冰川 累计 收获 了 42 根 有 效 
花 杆 ;2014 一 2015 年 对 Barpu 冰川 的 考察 计划 顺利 ， 
累计 收获 了 26 根 有 效 花 杆 ,2016 年 布设 的 花 杆 由 
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于 后 期 科 考 时 间 冲 突 , 未 能 完成 对 花 杆 的 后 期 再 观 
测 。 本 研究 中 的 花 杆 观测 信息 详 见 表 1、 表 2 和 表 3。 
同期 的 气温 数据 来 自 研 究 区 内 架设 的 自动 气象 
站 (automatic weather station,AWS)。 各 自动 气象 站 
均 在 距 下 热 面 2.0 mm 处 安装 有 温度 传感器 , 并 进行 
气温 数据 的 自动 观测 ,观测 间隔 为 30 min。 冰 川 区 
的 各 消融 点 的 气温 根据 冰川 上 的 气温 递减 率 插值 获 
得 。 插 值 的 基准 为 研究 区 内 各 冰川 获得 的 AWS 站 
点 数据 。 由 于 缺乏 2016 年 Sachen 冰川 消融 期 气温 
观测 数据 ,2016 年 度 的 气温 插值 基准 序列 由 上 一 年 
度 同 期 观测 序列 填补 。 温 度 递 减 率 利 用 站 点 气温 与 
其 对 应 的 海拔 高 度 , 采 用 线性 相关 方法 计算 获得 。 
其 中 Gharko 冰川 和 Barpu 冰川 的 温度 递减 率 是 根 
据 研 究 区 不 同 海拔 点 的 AWS 站 点 观测 资料 计算 求 
得 ,分 别 为 0. 92 % + (100m) !,0. 75 © 
(100m) ^^, Sachen 冰川 研究 区 内 只 架设 了 一 台 
AWS 气象 站 ,无 法 直接 计算 冰川 区 温度 递减 率 ,其 
温度 递减 率 是 利用 多 个 站 点 【包括 Sachen 冰 川 


R1 2014—2016 年 Gharko 冰川 花 杆 观测 点 数据 信息 
Tab.1 Information of the sighting rods for observing the Gharko Glacier from 2014 to 2016 


观测 年 度 海拔 /m 表 磺 厚度 /em 花 杆 数量 观测 时 段 人 月 -日 ) 有 效 观 测 次 数 
2014 3 346 ~3 502 7.8 ~30 ï 06 -06—07 -26 11 
2015 3 325 ~3 549 0.7 ~16.5 10 05 -30—08 -25 10 
2016 3 325 ~3 564 0~17 11 06 - 02 一 07 -26 11 
3 338 ~3 843 0-8 10 07 -23—10 -20 10 


#2 2014—2016 年 Sachen 冰川 花 杆 观测 点 数据 信息 
Tab.2 Information of the sighting rods for observing the Sachen Glacier from 2014 to 2016 


观测 年 度 海拔 /m 表 厂 厚 度 /em 花 杆 数量 观测 时 段 /( 月 - 日) ”有效 观测 次 数 
2014 3 587 ~3 861 1 ~20 9 06 - 16—07 - 19 9 
07 - 19—08 - 17 9 
2015 3 545 ~3 824 7 «19 11 06 - 15—07 -30 11 
07 -30 一 09 - 15 11 
2016 3 590 ~3 854 2 ~18.5 9 06 -30 一 08 - 11 9 
08 - 11—10 -28 8 
表 3 2014—2015 年 Barpu 冰川 花 杆 观测 点 数据 信息 


Tab.3 Information of the sighting rods for observing the Barpu Glacier from 2014 to 2015 


观测 年 度 海拔 /m 伦 杆 数量 表 太 厚度 /em 观测 时 段 人 月 -日 ) 有 效 观 测 次 数 
2014 2 939 ~3 002 13 2~10.5 07 - 11—08 -06 13 
2015 2 938 ~3 057 8 0 ~18 06 - 08 一 06 - 13 8 
2 ~23 06 - 13 一 08 - 16 4 
3 352 ~3 388 4 0~4.5 06 - 11 一 08 - 17 4 
3 559 1 1 10 -02 一 11 -07 1 
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表 4 冰川 观测 区 各 自动 气象 站 主要 观测 信息 
Tab.4 Main information of the automatic meteorological stations for observing glaciers 
冰川 序号 数据 时 段 海拔 温度 观测 高 度 m ”数据 采集 器 数据 采样 频率 /s 。 数据 记录 频率 /min 
Gharko AWSI 2014 -06 一 2016 -07 3 831 2 CR1000 10 30 
AWS2 2014 - 08 一 2017 -08 4 070 2 CR1000 10 30 
Sachen AWS 2014 -06 一 2016 -03 3 531 2 CR1000 10 30 
Barpu AWSI 2014 - 06 一 V2014 - 09 3 200 2 CR1000 10 30 
AWS2 2014 —06—2016 -08 2 784 2 CR1000 10 30 


AWS 气象 站 、 邻 近 Sachen 冰川 的 两 个 低 海 拔 气象 
站 Astore 气象 站 (2 167 m) 和 Rama 气象 站 (3 140 
m) ,其 中 Astore 气象 站 和 Rama 气象 站 距离 Sachen 
冰川 直线 距离 约 20 km ) 气温 数据 与 其 对 应 海拔 的 
线性 关系 ,通过 计算 斜率 获得 , 其 温度 递减 率 为 
0.65 C+ (100 m). 研究 区 内 自动 气象 站 布设 
信息 详 见 表 4。 

13 度 日 因子 计算 方法 


测 资料 时 发 现 :在 裸 冰 区 ,冰川 消融 量 与 正 积温 的 线 
性 关系 良好 ;但 是 在 表 矿 获 盖 区 ,其 线性 关系 不 甚 理 
想 。 本 研究 中 除 Barpu 冰川 外 ,Sachen 冰川 和 Ghar- 
ko 冰川 获得 的 冰川 消融 量 与 正 积温 之 间 的 线性 趋 
势 关 系 也 不 是 特别 理想 ,这 可 能 与 观测 资料 的 表 矶 
厚度 差异 分 布 有 关 ( 图 2)。 由 于 野外 考察 期 间 布 设 
的 花 杆 皆 位 于 冰川 的 消融 区 ,而 消融 区 覆盖 着 大 量 


度 日 模型 是 基于 雪 冰 消 柄 量 与 气温 关系 的 统计 500 MET 
模型 ,通过 建立 度 日 因子 与 正 积温 之 间 的 数学 函数 -| m ee 
关系 推 求 冰川 消融 量 ""-"”。 根 据 定义 , 度 日 因子 i oso | ——-— 
的 基本 表现 形式 如 下 : E 2000 上 m uu 
M = PDD x DDF (1) x 1 000 1 2 
式 中 :MM 为 某 一 时 段 冰 川 或 雪 的 消融 水 当量 ,mm w. 500 1000 —— 1500 
e. ;DDF 为 冰川 或 雪 的 度 日 因子 ,mm «d^! Cy 
. 6000 - (CDb)Gharko 冰 川 
PDD 为 某 一 时 段 内 的 正 积温 ,其 一 般 由 下 式 获取 : S " 
Vot 
PDD = 5H, .了 (2) : 4000 f Pu Sec ey 
、 a4 Pl d ?3000 | 2-77 * y2.8262x+723 
SUP HH, AZE, T, 20 BE,H, 21.0, 34 T, 0 & 200 Qe 。。 R?=0.251 4 
时 ,及 20.0; T, 为 某 天 (1 的 日 平均 气温 (% ) 。 = 
a 1 
2 结 果 5i tit p 700 E i S 100 1300 1500 
Me 3500 - (@Sachen 冰 川 
气温 是 反映 冰川 热量 综合 状况 的 良好 指标 , 且 300 | women t 
容易 通过 气象 站 数据 与 温度 递减 率 插值 获得 国内 外 ve e ete co 
外 众多 研究 结果 表明 ,冰川 消融 与 正 积温 之 间 存 在 ah 150| op A 。 
良好 的 线性 关系 。 度 日 因子 正 是 在 此 基础 上 建立 。 000 | UM 
的 。Braithwaite 在 对 格陵兰 冰 盖 研究 中 发 现 , 冰 a I [Oo te nM 
川 消融 量 与 正 积温 的 相关 系数 达 0.9 以 上 ;在 国内 ， a å 


冰川 工作 者 们 对 念 青 唐 古 拉 山 扎 当 冰川 ,祁连山 
“七 一 ”冰川 和 唐古拉 山 冬 克 玛 底 冰 川 ”的 研究 
中 ,也 发 现 了 冰川 消融 与 气温 之 间 的 良好 趋势 关系 。 
然而 ,并 不 是 所 有 观测 研究 结果 都 符合 上 述 现象 。 
张 勇 等 ”分析 在 天 山南 坡 科 其 卡尔 巴 契 冰川 的 观 


图 2 Sachen 冰川 .Gharko 冰川 和 Barpu 冰川 
消融 量 与 正 积温 的 关系 
Fig.2 Relationship between the ablation amount and the 
cumulative positive temperature of the Sachen, Gharko 


and Barpu glaciers 


Pa 
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的 差异 。 其 中 Barpu 冰川 观测 点 的 表 磺 厚度 以 小 
F10 cm 厚度 点 为 主 ( 图 3)。 与 Barpu 冰川 相 比 ， 
Sachen 冰川 和 Gharko 冰川 的 观测 点 表 厂 厚度 总 体 
WA , 表 矿 厚度 介 于 10 ~ 20 em 之 间 的 观测 结果 较 
多 。 
2.1 度 日 因子 随时 间 变 化 特征 

利用 花 杆 记录 的 冰川 消融 数据 各 自动 气象 站 
逐日 气温 数据 和 温度 递减 率 推 求 各 花 杆 监测 点 的 冰 
川 度 日 因子 。 需 要 特别 说 明 的 是 ,由 于 缺乏 2016 年 
Sachen 冰川 消融 期 AWS 气温 观测 数据 , 故 2016 年 
Sachen 冰川 各 花 杆 观测 点 的 消融 期 温度 数据 是 基 
于 2015 年 AWS 记录 的 同期 气温 数据 插值 获得 的 。 
根据 式 (1) 和 式 (2) 计算 并 汇总 度 日 因子 数据 结 

( 表 5)。 计 算 结果 显示 :Sachen 冰川 上 的 各 观测 时 
段 布设 的 花 杆 点 度 日 因子 均值 介 于 2. 40 ~ 3. 16 
mm + d`! .C7 ,总 体 均值 2.83 mm d - C753 
Gharko 冰川 上 的 各 观测 时 段 布设 的 花 杆 点 度 日 
子 均值 介 于 2.97 ~5.41 mm. d^! - C 7! , AGES E 
3.74 mm . d! + °C ^! ; Barpu 冰川 上 的 各 观测 时 段 
布设 的 花 杆 点 度 日 因子 均值 介 于 2.75 ~7.09 mm * 
d' C, RUA 3.91 mm d + "C7, BA 
因子 在 Gharko 冰川 和 Sachen 冰川 变化 幅度 较 大 ， 
而 Sachen 冰川 度 日 因子 较为 稳定 。 受 冰 面 气象 要 
素 、 水 热 条 件 以 及 地 形 因 素 的 影响 ,冰川 上 的 能 量 交 
换 过 程 会 随 着 时 间 变 化 而 变化 ,因此 ,各 冰川 研究 区 
度 日 因子 在 时 间 尺 度 上 存在 差异 ,具体 表现 为 :同一 
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图 3 Sachen 冰川 .Gharko 冰川 和 Barpu 冰川 度 日 因子 
PERDUT REZE 
Fig.3 Variation of degree-day factors of the Sachen , Gharko 


and Barpu glaciers with the thickness of surface moraine 


#5 Sachen 冰川 .Gharko 冰川 和 Barpu 冰川 度 日 因子 计算 结果 
Tab.5 Calculated results of degree-day factors of the Sachen,Gharko and Barpu glaciers 


冰川 年 份 观测 时 段 /月 -日 各 花 杆 观测 点 度 日 因子 计算 结 度 日 因子 均值 年 平均 值 总 体 均值 
Sachen 2014 06 - 16—07 - 19 1.14 ~4.70 3.16 2.78 2.83 
07 - 19 一 08 - 17 1.09 ~5.73 2.40 
2015 06 - 15—07 -30 0.89 ~4.22 2.98 2.81 
07 -30—09 -15 1.77 ~3.89 2.88 
2016 06 -30 一 08 - 11 2.17 ~5.28 3:13 2.90 
08 - 11—10 -28 1.84 ~3.45 2.63 
Gharko 2014 06 -03 一 07 -26 1.02 ~4.49 2:97 2.97 3.74 
2015 05 -30—08 -25 2.74 ~4.16 3.51 3.51 
2016 06 -02—07 - 26 1.54 «5.63 3.22 4.26 
07 -24—10 -20 3.01 «6.81 5.41 
Barpu 2014 07 - 11—08 - 06 2.09 ~5.53 3.45 3.45 3.91 
2015 06 -08—06 - 13 1.34 ~6.22 4.45 4.27 
06 - 11—08 - 17 4.48 ~4. 84 4.73 
06 - 13—08 - 16 2.08 ~3.43 2.75 
10 -02—11 -07 7.09 7.09 
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图 4 2015 年 Sachen 冰川 布设 的 11 根 花 杆 不 同 观测 时 段 
冰川 度 日 因子 变化 


Fig.4 Variation of degree-day factors of the Sachen glacier 
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度 日 因子 /mm * d^ * 


9 10 11 


observed from 11 sighting rods during different observational 


periods in 2015 


花 杆 观测 点 的 度 日 因子 在 不 同时 间 上 的 计算 结果 不 
尽 相 同 ,例如 2015 年 布设 于 Sachen 冰川 表面 的 11 
根 花 杆 在 6 月 15 日 至 7 月 30 日 与 7 月 30 日 至 9 
月 15 日 的 2 个 连续 观测 时 段 内 计算 的 度 日 因子 值 
都 存在 不 同 程度 的 差异 (图 4) ,这 说 明 冰 川上 单 点 
的 度 日 因子 会 随时 间 的 变化 而 变化 ; 随 着 观测 时 段 
的 变化 ,基于 各 观测 时 段 计 算 的 度 日 因子 均值 出 现 
明显 的 波动 ;各 冰川 度 日 因子 的 年 平均 值 波 动 幅 度 
较 小 , 旦 逐年 增加 的 趋势 。 
2.2 度 日 因子 的 空间 分 布 特征 

由 于 度 日 因子 会 随 着 时 间 的 变化 而 变化 , 花 杆 
观测 点 在 不 同 观测 时 段 下 的 度 日 因子 计算 结果 反映 
了 单 点 度 日 因子 在 时 段 内 的 均值 水 平 , 因 而 较 长 时 
段 下 计算 的 度 日 因子 值 更 适合 反映 花 杆 观测 点 度 日 
因子 的 实际 “平均 状态 ”"。 基 于 上 述 认识 ,笔者 得 到 
了 各 冰川 研究 区 花 杆 观测 点 度 日 因子 与 海拔 的 对 应 
关系 ,如 图 5 所 示 (Sachen 冰川 各 观测 花 杆 存在 两 
个 连续 时 段 的 度 日 因子 计算 结果 ,由 于 前 后 观测 时 
段 时 长 较 接近 , 故 取 其 均值 作为 各 花 杆 点 度 日 因子 
值 ;Barpu 冰川 2015 年 布设 的 首 批 8 根 花 杆 只 有 4 
根 获 得 了 较 长 时 段 的 观测 记录 ,因此 选取 这 4 次 观 
测 结果 与 其 他 花 杆 点 度 日 因子 作为 本 节 研 究 数 
据 )。 由 图 5 可 以 看 出 , 度 日 因子 在 空间 上 存在 明 
显 的 差异 ,整体 上 度 日 因子 存在 随 海拔 高 度 增加 而 
递增 的 趋势 ,这 一 变化 趋势 与 学 者 们 对 青藏 高 原 地 
KREDKI BOUES A a BO 
分 析 了 中 国 西部 数 十 条 冰川 的 空间 变化 特征 ,发现 
较 大 的 度 日 因子 通常 出 现在 海拔 较 高 的 地 区 ,尤其 
是 冰川 物质 平衡 线 附 近 。 这 主要 是 冰川 物质 平衡 线 
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Fig.5 Variation of degree-day factors of the Sachen, 


Gharko and Barpu glaciers with the altitude 


附近 太阳 辐射 强烈 ,导致 冰川 消融 强烈 造成 
的 。 此 外 ,还 有 研究 表明 冰川 物质 平衡 线 附 近 
存在 所 谓 的 低温 效应 ,其 气温 较 低 ,因而 正 积温 量 值 
BENNO? 。 本 研究 中 涉及 的 各 花 杆 观测 点 布设 于 冰 
川 消融 区 , 随 着 海拔 的 增加 , 愈 靠近 冰川 物质 平衡 
线 , 度 日 因子 愈 呈 递增 趋势 ; 受 区 域 差 异 影响 ,不 同 
冰川 的 度 日 因子 空间 递增 率 不 同 。 基 于 各 观测 冰川 
所 有 观测 点 的 度 日 因子 值 ,笔者 建立 了 Sachen 冰 
川 、Gharko 冰川 和 Barpu 冰川 度 日 因子 与 海拔 间 的 
拟 合 函数 关系 ,分 别 为 : 

Sachen 冰川 : 

DDF = 0.003 7H - 10.962,( 尼 = 0.324 7) 

Gharko 冰川 : 

DDF =0.007 4H - 21.92,(R? = 0.499 9) 


Barpu 冰川 : 
DDF = 0.004 1H - 8.914 5,(R? = 0.418 7) 
式 中 :DDF 为 冰川 度 日 因子 (mm + d^! - C7!) 5H 
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度 递减 率 较 大 , 正 积 温 量 值 变化 偏 小 造成 面向 阳 坡 
的 冰川 度 日 因子 随 海拔 递增 率 较 大 。 

对 比 青 藏 高 原 西 部 冰川 度 日 因子 的 研究 成 果 
( 表 6) ,发 现 Sachen ZKJI| , Gharko 冰川 和 Barpu vk 


为 海拔 高 度 (m) 。 可 由 上 和 式 获 得 各 海拔 高 度 处 的 度 
日 因子 值 ,从 而 开展 区 域 分 布 式 冰川 物质 平衡 和 融 
水 径流 模拟 研究 。 

对 比 冰川 朝向 会 发 现 , 朝 向 为 SW 的 Gharko vk 
川 度 日 因子 随 海拔 的 递增 率 明 显 大 于 朝向 同 为 NE 
的 Barpu 冰川 和 Sachen 冰川 ,在 高 海拔 区 域 Gharko 
冰川 度 日 因子 偏 大 ,这 与 韩 海 东 等 在 科 其 喀 儿 冰 
川 表 厂 区 研究 的 不 同 坡 向 冰 岸 度 日 因子 变化 规律 一 
致 ,说 明度 日 因子 变化 受 冰 川 朝向 影响 ,面向 阳 坡 的 
Gharko 冰川 在 消融 期 接收 的 太阳 辐射 较 大 ,冰川 消 
融 更 强烈 ,同时 较 高 海拔 地 区 温度 较 低 且 该 冰川 温 


川 度 日 因子 均值 与 青藏 高 原 西部 冰川 度 日 因子 均值 
较 接 近 。 张 勇 等 2 根据 冰川 度 日 因子 的 均值 ,将 中 
国 西部 冰川 区 分 为 3 个 区 域 : 度 日 因子 高 值 区 
(29.0mm * d^ * C7), 度 日 因子 中 值 区 
(6.0~9.0 mm * d! + ^C 7) 和 度 日 因子 低 值 区 
(<6.0 mm-d' C7) 。 本 文 研究 的 Sachen, 
Gharko 和 Barpu 冰川 所 在 的 喀 喇 昆仑 一 西 喜马拉雅 
山脉 区 属于 其 划分 的 度 日 因子 低 值 区 ,各 冰川 研究 
区 度 日 因子 均值 计算 结果 丝 小 于 6.0 mm - 
dC ,与 其 研究 成 果 吻 合 。 


表 6 印度 河上 游 流域 及 其 周边 地 区 冰川 度 日 因子 对 比 


Tab.6 Compared results of the degree-day factors of glaciers in the upper Indus and its surrounding areas 


冰川 度 日 因子 /(mm: d^! - € 7!) 观测 时 段 海拔 /m 参考 文献 

巴 托 拉 冰 川 3.4 1975 -06 -08 2 780 29 
特 拉 木 坎 力 冰川 5.9 1987 -06 -25 一 09 -07 4 630 29 
6.4 1987 -06 -24 一 09 -07 4 650 29 
切 尔 干 布 拉 克 冰川 2.6 1960 -06 -06—07 -30 4 750 29 
洋 布 拉克 冰川 4.3 1987 -07 -01—07 -05 4 800 29 
Mg" vl 3.0 1989 -05 -07—09 - 07 4 840 29 
Sek 4.7 1981 - 06 - 15—06 - 28 4 616 29 
3.6 1987 -06 - 14—06 -27 4 900 29 

Sachen 冰川 2.83 本 研究 野外 考察 时 段 3 526 ~3 861 本 研究 

Barpu 冰川 3.91 本 研究 野外 考察 时 段 2 938 ~3 576 本 研究 

Gharko 冰川 3.74 本 研究 野外 考察 时 段 3 325 -3 843 本 研究 


2.3 ” 表 磺 对 度 日 因子 影响 分 析 
度 日 因子 是 冰川 表面 能 量 收文 状况 的 宏观 表 


冰川 消融 的 主要 能 源 是 太阳 辐射 。 在 热量 
收入 中 , 耗 于 冰 面 消融 的 热量 与 冰 面 性 质 有 关 ,对 裸 


现 。 根 据 不 同 研究 者 对 度 日 因子 的 研究 结果 , 度 日 
PIF SE yk ER LE wg 0006280090 Be qr Jer 
FEO? p BRUT 地 形 ' 等 自然 因素 的 变化 都 会 
对 冰川 消融 产生 一 定 的 影响 ,进而 影响 度 日 因子 的 
变化 。 此 外 , 花 杆 人 工 测 量 以 及 基于 气象 站 插值 生 
成 的 温度 数据 带 来 的 不 确定 性 也 不 容 忽视 。 换 
而 言 之 ,冰川 消融 点 的 度 日 因子 计算 结果 受 自 然 和 
人 为 因素 综合 影响 。 由 于 笔者 在 野外 科 考 中 仅仅 获 
得 了 各 花 杆 点 的 表 磺 厚度 变化 数据 , 除 表 磺 因素 外 
的 其 他 影响 因子 在 本 文中 皆 无 法 开展 深入 探讨 , 故 
聚焦 于 表 古 变化 对 度 日 因子 的 影响 研究 。 


露 冰 面 而 言 ,90% 以 上 的 热量 供 冰 面 消 融 “”。 表 三 
的 覆盖 作用 改变 了 冰川 与 大 气 之 间 的 热量 传输 过 
程 。 与 裸 冰 区 相 比 , 表 厂 的 覆盖 作用 削弱 了 地 面 反 
射 率 ,促进 了 地 面 辐射 的 吸收 。 而 地 面 辐射 只 能 先 
加 热 表 厂 履 盖 层 ,然后 再 将 热量 传导 至 冰川 表面 。 
因此 , 表 矿 的 存在 会 影响 到 达 并 用 于 冰 面 消融 的 耗 
热量 。 相 关 研 究 表明 "” : 表 厂 厚度 分 布 的 异 质 性 会 
加 速 或 抑制 冰川 消融 这 一 过 程 ,从 而 深刻 影响 冰川 
消融 过 程 , 当 表 矿 层 厚度 小 于 某 一 临界 厚度 时 ,其 下 
覆 冰 川 冰 消 融 大 于 裸 冰 消融 , 即 表 厂 覆盖 加 速 冰 川 
消融 ; 随 着 表 矿 厚度 的 增 大 (大 于 临界 值 ) , 表 矿 的 
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阻 热 作 用 抑制 了 冰川 消融 。 度 日 因子 是 观测 时 段 内 
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作用 促进 了 冰川 的 消融 , 表 奔 覆盖 区 冰川 度 日 


的 冰川 消融 量 与 同期 正 积 温 的 比值 ,不 同 厚度 的 表 
奢 物 对 冰川 消融 的 差异 影响 势必 会 导致 度 日 因子 的 
变化 。 

由 于 不 同 观测 时 段 获 得 的 花 杆 观测 消融 量 对 
应 的 正 积 温 量 差异 较 大 ,基于 正 积温 计算 的 观测 
期 日 平均 气温 差异 也 显著 ,此 外 ,在 野外 观测 中 布 
设 的 裸 冰 面 花 杆 点 较 少 ,为 了 有 效 探讨 裸 冰 及 不 
同 表 奢 厚度 影响 下 的 冰川 消融 量 与 度 日 因子 的 关 
系 , 笔 者 选择 Barpu 冰川 2015 年 6 月 8 一 13 日 观 
测 时 段 内 布设 8 根 花 杆 和 Gharko 冰川 2016 年 6 
月 2 日 至 7 月 26 日 观测 时 段 内 布设 11 根 花 杆 观 
测 资料 做 初步 探讨 (图 6a 和 图 6b) 。 从 图 6a 和 图 
6b 可 以 看 出 , 表 矿 厚度 低 于 2 cm 时 , 表 磺 的 覆盖 
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表 研 厚度 /cm 
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子 大 于 裸露 区 冰川 的 度 日 因子 。 在 Sachen 冰川 由 
于 野外 花 杆 观测 工作 中 没有 布设 裸 冰 面 下 的 花 杆 
点 , 故 无 法 直接 判明 Sachen 冰川 上 的 表 磺 对 度 日 
因子 影响 的 临界 厚度 值 ,本 文 不 对 Sachen 冰川 的 
表 夸 影响 做 深入 探究 。 需 要 特别 说 明 的 是 :冰川 
度 日 因子 对 应 的 表 奢 厚度 是 指 观测 时 段 的 前 后 2 
次 观测 记录 中 获得 的 表 奢 厚度 的 平均 值 。 然 而 ， 
大 多 数 花 杆 观 测 工 作 的 间隔 时 间 较 长 , 受 地 表 外 
力 和 冰川 运动 影响 , 冰 面 表 磺 物 可 能 在 冰川 表面 
发 生 移动 , 表 磺 厚度 变化 存在 一 定 的 不 确定 性 , 故 
本 研究 中 使 用 的 平均 表 奔 厚度 的 可 靠 度 有 待 进 一 
步 商 榨 。 基 于 以 上 认识 ,本 文 不 对 表 磺 物 对 冰川 
消融 影响 的 临界 厚度 值 做 出 结论 。 
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注 :(a) Barpu 冰川 2015 年 6 月 8 一 13 日 观测 时 段 布设 的 8 根 花 杆 观测 点 ,(b) Gharko 冰川 2016 年 6 月 2 日 至 7 月 26 日 
观测 时 段 布设 的 11 根 花 杆 观测 点 。 
图 6 不 同 表 硬 厚度 影响 下 永川 日 平均 消融 速率 日 平均 气温 和 度 日 因子 变化 


Fig.6 Variation of daily glacier melt rate, daily average temperature and degree-day factors with different thickness 


of surface moraine 


各 冰川 研究 区 表 研 广 布 , 厚 度 不 一 ,主要 集中 在 
消融 区 。 其 中 Sachen 冰川 表 研 覆 盖 率 最 高 , 除 少量 
海拔 区 域外 ,多 被 表 厂 覆盖 ;Charko 冰川 次 之 ,Barpu 
冰川 最 少 。 本 文 涉及 的 各 冰川 研究 区 消融 点 度 日 
子 与 冰川 表面 表 厂 厚度 的 对 应 关系 (图 3)。 从 图 3 
可 以 看 出 ,各 冰川 研究 区 度 日 因子 与 表 研 厚度 之 间 


表现 出 比较 一 致 的 梯度 变化 ,整体 上 度 日 因子 随 表 
厂 厚 度 的 增加 而 递减 , 仅 表 厂 厚度 小 于 2 cm 的 部 分 
冰川 区 表现 出 相反 的 梯度 变化 (由 于 较 小 的 表 厂 厚 
度 在 一 定 程度 上 促进 了 冰川 消融 ,致使 度 日 因子 增 
加 )。 基 于 各 观测 冰川 所 有 观测 点 的 度 日 因子 值 ， 
初步 建立 了 Sachen 冰川 .Gharko 冰川 和 Barpu 冰川 
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度 日 因子 与 表 厂 厚度 间 的 整体 线性 拟 合 函 数 关 系 ， 了 一 定 的 影响 ,朝向 为 SW 的 Charko 冰川 度 日 因子 

分 别 为 : 随 海拔 递增 率 明 显 大 于 朝向 同 为 NE 的 Barpu 冰川 
Sachen 冰川 : 和 Sachen 冰川 。 

DDF = - 0.110 5X + 4.009 9,(R® = 0.286 5) (4) 度 日 因子 受 表 磺 变 化 影响 显 着 。 整 体 上 上 度 
Charko 冰川: 日 因子 随 表 矿 厚 度 的 增加 而 递减 ,在 表 厂 厚度 小 于 

Mcr 2 cm 时 , 表 厂 的 覆盖 作用 促进 了 冰川 的 消融 , 表 矿 

i ` 覆盖 区 冰川 度 日 因子 大 于 裸露 区 冰川 的 度 日 因子 。 

Barpu 冰川 : 随 着 表 厂 厚 度 的 增加 , 表 厂 物 逐 渐 前 弱 了 冰 面 消融 

DDF = - 0.182X + 4.981 2,(R - 0.590 8) 的 耗 热量 ,冰川 度 日 因子 总 体 上 呈 递 减 趋 热 。 

式 中 :DDF 为 冰川 度 日 因子 (mm d + Co!) XH 


表 厂 厚度 (cm) 。 可 由 上 式 估 算 不 同 表 厂 厚 度 下 的 
度 日 因子 值 。 
3 结语 


本 研究 利用 印度 河上 游 流 域 Sachen 冰川 、 
Gharko 冰川 和 Barpu 冰川 实测 花 杆 和 自动 气象 站 记 
录 气 温 数 据 资 料 估算 了 各 冰川 研究 区 观测 花 杆 点 的 
度 日 因子 值 ,分 析 了 冰川 度 日 因子 的 在 时 间 和 空间 
上 的 变化 特征 ,初步 探究 了 度 日 因子 的 影响 因素 ,可 
以 得 到 以 下 结论 : 

(1) 计算 结果 显示 :2014 一 2016 年 Sachen 冰川 
度 日 因子 为 0.89 ~5.73 mm - d^! - 7C ,各 观测 时 
段 度 日 因子 均值 介 于 2.40 ~3.16 mm d'! - C75, 
总 体 均 值 为 2.83 mm . d^! + C7!; Gharko 冰川 
2014—2016 年 度 日 因子 为 1.02 ~ 6. 81 mm - 

A acu dad g- 2:97 
5.4] mm: d! - «^! , jJ 5 (E 3. 74 mm - 
d^! + °C 7 ;2014—2015 年 Barpu 冰川 上 观测 获得 的 
度 日 因子 介 于 1.34 ~7.09 mm . d^! - Co! AI 
时 段 度 日 因子 均值 为 2.75 ~7.09 mm - d - Co, 
总 体 均 值 为 3.91 mm +d) + C^, 

(2) 在 时 间 分 布 上 , 度 日 因子 会 随 着 时 间 变 化 
而 改变 。 不 同 观测 时 段 计 算 的 度 日 因子 值 存在 明显 
差异 ,即使 是 同一 观测 点 的 度 日 因子 在 不 同 观 测 时 
段 的 计算 结果 也 不 同 。 但 是 度 日 因子 年 平均 值 年 际 
变化 幅度 不 大 , 呈 逐 年 递增 的 趋势 。 

(3) 在 空间 分 布 上 ,冰川 度 日 因子 空间 变化 主 
要 体现 为 随 海拔 与 地 形 因素 的 变化 而 变化 。 主 要 表 
现 为 :Sachen 冰川 .Gharko 冰川 和 Barpu 冰川 度 日 
子 随 海拔 的 上 升 而 增 大 , 随 海拔 的 递增 率 分 别 为 
0.003 7 mm * d^! - C ~! 0.0074 mm- d^! - C7" 
0.004 1 mm * d^' + ^C ;冰川 朝向 对 度 日 因子 产生 
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Abstract: ^ Currently ,the degree-day model based on the sum of daily mean temperature is widely used in glacier 
ablation research. The degree-day factors for ice ablation of the selected Sachen, Gharko and Barpu glaciers were 
calculated based on the data of mass balance and meteorology from 2014 to 2016 , and their characteristics including 
the spatiotemporal variation and affecting factors were analyzed. The results showed that the average degree-day fac- 
tors of the Sachen , Gharko and Barpu glaciers were 2. 83,3. 74 and 3.91 mm * d! - C ^! respectively. The degree- 
day factors of the glaciers increased with the increase of elevation , and the increase rates of degree-day factors of the 
Sachen , Gharko and Barpu glaciers were 0. 003 7,0. 007 4 and 0. 004 1 mm * d! - C ^! respectively. The degree- 
day factor was not a constant, it changed with the time and was affected significantly by surface moraine , and it was 
holistically decreased with the increase of surface moraine thickness. The coverage of surface moraine thinner than 2 
cm could promote the glacier ablation, and the degree-day factor in the glacier area covered by surface moraine was 
higher than that without surface moraine. Glacier aspect affected the degree-day factor to a certain extent, and the 
degree-day factor on the sunny slope was higher than that on the shady slope. 
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